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Abstract 

Evolution of an anisotropic universe described by a Bianchi type I (BI) model 

\ in presence of nonlinear spinor field has been studied by us recently in a 

I series of papers. On offer the Bianchi models, those are both inhomogeneous 

' and anisotropic. Within the scope of Bianchi type VI (BVI) model the self- 
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consistent system of nonlinear spinor and gravitational fields are considered. 
The role of inhomogeneity in the evolution of spinor and gravitational field is 
studied. 
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1. INTRODUCTION 



Several authors studied the nonlinear spinor fields (NLSF) since Ivanenko |T]-0] showed 
that a relativistic theory imposes in some cases a fourth order self-coupling. Nonlinear spinor 
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field in an anisotropic universe, namely in a Bianchi-type I universe (BI) is studied by us in 
a series of papers p]-|TT| . In these papers we considered the nonlinear spinor field, as well as 
a system of interacting spinor and scalar fields. Beside the spinor fields, we also study the 
role of a A term in the evolution of the universe. For the details of a NLSF in a BI universe 



one can consult [12]. In the light of results obtained in the previously mentioned papers 
the study of NLSF in other anisotropic universes presents great interest. In this report we 
consider the self-consistent system of nonlinear spinor and BVI gravitational fields. The 
results have been compared with those obtained for the BI universe. 

2. FUNDAMENTAL EQUATIONS AND GENERAL SOLUTIONS 

We choose the Lagrangian for the self-consistent system of spinor and gravitational fields 
in the form 



R i 
2k + 2 



— Mipip + Ln (2.1) 



with R being the scalar curvature and k being the Einstein's gravitational constant. The 
nonlinear term L n describes the self-interaction of spinor field and can be presented as some 
arbitrary functions of invariants generated from the real bilinear forms of spinor field having 
the form: 

S = Tpip (scalar), (2.2a) 
P = i'tpj ip (pseudoscalar), (2.2b) 
= (frfi{>) (vector), (2.2c) 
A 1 * = (?/>7VV) (pseudovector), (2. 2d) 

QT = (ipa^-ip) (antisymmetric tensor), (2.2e) 

where a^ u = (i/2)[-f' J ''-f u — 7^7^]. Invariants, corresponding to the bilinear forms, are 

/ = S 2 , (2.3a) 
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J = P 2 , (2.3b) 

I v =v tt v»= (^ 7 M V0 VWVO, (2-3c) 

I A = ApA»= (V^VVO M^VVV). (2.3d) 

Iq = Qnv = $<tVi/>) g^g v 0a al3 ijj). (2.3e) 

According to the Pauli-Fierz theorem |]T3| among the five invariants only / and J are 

independent as all other can be expressed by them: Iy = —Ia = I + J and Iq = I — J. 



Therefore, we choose the nonlinear term Ln = F(I, J), thus claiming that it describes the 
nonlinearity in the most general of its form. 

We choose the anisotropic inhomogeneous universe given by a Bianchi-type VI (BVI) 
metric: 

ds * = dt 2 _ a 2 e -2 mz dx 2 _ b 2 e 2nz dy 2 _ Q 2 ^ q q 

with a, b, c being the functions of time only. Here m, n are some arbitrary constants and 
the velocity of light is taken to be unity. Note that the suitable choice of m, n as well 
as the metric functions a, b, c in the BVI given by (|2.4j ) evokes the following Bianchi-type 
universes. Thus 

• for m = n the BVI metric transforms to a Bianchi-type V (BV) one, i.e., 
m = n, BVI =^ BV G open FRW; 

• for n = the BVI metric transforms to a Bianchi-type III (Bill) one, i.e., 
n = 0, BVI Bill; 

• for m = n = the BVI metric transforms to a Bianchi-type I (BI) one, i.e., 
m = n = 0, BVI BI; 

• for m = n = and equal scale factor in all three direction the BVI metric transforms 
to a Friedmann- Robertson- Walker (FRW) universe, i.e., 

m = n = and a = b = c, BVI FRW. 
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Variation of the Lagrangian ( |2.1|) with respect to field functions ip{ip) gives the nonlinear 
spinor field equations: 

n^V^V — Mip + Vifj + iG^ip = (2.5a) 

Wurf'f + — — i^7 5 = (2.5b) 

where V = 2SF, and Q = 2PFj. 

Varying ( |2.1| ) with respect to metric function (g^) we find the Einstein equations for 
a, b, c which for the metric ( |2.4| ) read 

s + - + r - ^ = P-fc) 

c oc <r 
c a ca m 2 2 

- + - + =- = -kT 2 , (2.6b) 

c a ca cr 

a 6 d 6 mn ~ . 

- + 7 + -T + — = -^|, 2.6c 
a b ab cr 

ab be ca m 2 - mn + n 2 , n 

~7 + I~ + 5 = - kT 0) 2 - 6d 

ab be ca cr 

m- — n- — (m — n)- = — kT2. (2.6e) 

a b y J c 6 K ' 

Here overdots denote differentiation with respect to time (t) and is the energy-momentum 
tensor of the material field and has the form 

T£ = \g pv (fcV^ + i'luV^ - V^7^ - V v $lJ) ~ ^L sp . (2.7) 

Here L sp is the spinor field Lagrangian, which on account of spinor field equations ( |2.5| ) takes 
the form: 

L sp = -VS -QV + F. (2.8) 
In the expressions above V„ denotes the covariant derivative of spinor, having the form ||14|: 



V M V = - (2.9) 
where T M (x) are spinor affine connection matrices defined by the equality 

r M (x) = (1/4)^(^(9^ - ryyV- 
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For the metric element (|2.4| ) it gives 

r = o, 

ri = \{ai V - m a -ff]e— 

— C 

T 3 = Uff 



It is easy to show that 



where we define 



, m It n m — n , 



(2.10) 



The Dirac matrices 7 M (x) of curved space-time are connected with those of Minkowski one 
as follows: 

7° =7°, 7 1 = 7 1 e m 7a, 7 2 = f/6e" 2 , 7 3 = f/c 

with 
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v°- J / 



( n .A 



a 1 
-a 1 



, 7 5 = 7 5 



—J 



where a,- are the Pauli matrices: 



a 



f° 











fl ] 














v 1 






^ J 







Note that the 7 and the a matrices obey the following properties: 

7¥ + f¥ = 2^', i,j = 0, 1, 2, 3 
ff+ff = 0, (7 5 ) 2 = /, z = 0,l,2,3 
a 3 a k = 5 jk + is jk icr l , j, fc,/ = l,2,3 

where 77^ = {1, —1, —1, —1} is the diagonal matrix, 8jk is the Kronekar symbol and Ejki is 
the totally antisymmetric matrix with £123 = +1. Let us consider the spinors to be functions 
of t and z only, such that 



Bijan Saha 

ip(t, z) = v(t)e ikz } $(t, z) = v{t)e~ ikz 



Inserting (|2.11| ) into (|2.5|) for the nonlinear spinor field we find 
f \ / m — n k 



- i- I7 v + z$f + £7 v = 
2c c> 

in — n k\ q _ „ . 

h i- t>7 - i®v - Qvy = 0. 

2c c 7 



(2.11) 

(2.12a) 
(2.12b) 



Here we define $ = M — T>. Introduce a new function Uj(t) = y/rvj(t), for the components 
of the NLSF from ( |2.12 ) one obtains 



111 + i®v>i — 

11 2 + i&u 2 + 
u 3 - i§u 3 - 

U4 — i§U4 + 



m — 


n 


2c 




m — 


n 


2c 




m — 


n 


2c 




m — 


n 


2c 



k 1 

- i- + Q u 3 = 0, 

c J 

k i 

— i Q m 4 = 0, 

c J 



-%--g 

c 

k „ 
-i- + Q 

c 



ui = 0, 
u 2 = 0. 



(2.13a) 
(2.13b) 
(2.13c) 
(2.13d) 



Using the spinor field equations Q2.5|) and ( [2.12|) it can be shown that the bilinear spinor 
forms, defined by ( p. 2D , i.e., 

S = tpip = vv, P = i^j ijj = iv^v, A = ^7 5 7°'?/> = v*/ *f°v , 

A 3 = ^7 5 7 3 '0 = vrfrfv, V° = ■07°'0 = ^7%, V 3 = ■^7 3 V ; = v*f 3 v, 

Q 30 = i^rfrfij) = w7 3 7%, Q 21 = ^7°7 3 7 5 '?/' = iip^ 2 ^ 1 ^ = iv^^v, 

obeying the following system of equations: 

0, 



k 



S — 2—Q — 2QA Q 
P - - 2$A° Q = 0, 

A° —A 3 + 2$P + 2GS = 0, 

• 3 m — n 
A - A = 0, 



1>° - 



c 

m — n. 



-V 3 = 0, 



U 3 
1 o 



m — n 



V ° + 2$Q 30 - 2QQf = 



k 



+ 2-S - 2$V 6 = 0, 
Ql 1 + 2^P + 2QV 3 = 0, 



(2.14a) 
(2.14b) 
(2.14c) 
(2.14d) 
(2.14e) 
(2.14f) 
(2.14g) 
(2.14h) 
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where we denote F = tF. Combining these equations together and taking the first integral 
one gets 



(So) 2 + (Po) 2 + (Ay - (AlY + (V o y + (VfY + (Qt") 2 + (QfY = C = Const (2.15a) 

Before dealing with the Einstein equations (|2.6|) let us go back to ( |2.13|) . From the first 
and the third equations of the system ( |2.13| ) one finds 



«ia = (Q- Q)u 2 13 - 2i$ Ul3 + (g + Q), 



(2.16) 



where, we denote Wi 3 = ui/u 3 and Q — [m — n — 2ik]/2c. The equation ( |2 . 1 6[ ) is a Riccati 
T5| with variable coefficients. Further, setting w 13 = w 13 exp[— 2i J $>(t)dt], from ( J2.16 ) 



one 



one obtains 



*>13 



(Q - Q)v 2 13 e- 2 *f •<*>* + (Q + Q)e 2l I 



(2.17) 



The general solution of ( [2. 17| ) can be written as 



Vl3 



J(G -Q)e~ 2i I^ dt dt + C(t) 



(2.18) 



with the integration constant be defined from 



r i 

/ J(G -Q)e~ 2i $^ t)dt dt + C(t) dC = J(G + Q)e 2i $^ t)dt dt. (2.19) 

Thus given the nonlinear term in the Lagrangian and solution of the Einstein equations, one 
finds the relation between u\ and u% (u2 and as well), hence the components of the spinor 
field. 



Now we study the Einstein equation ( |2.6| ). In doing so, we write the components of the 
energy-momentum tensor, which in our case read 



T ° = mS - F + -V 3 , 



Tl = T 2 = VS + QP — F. 

k 

c 



Tn = VS + QP — F — -V 3 



T 3 ° = -kV°. 



(2.20a) 
(2.20b) 
(2.20c) 
(2.20d) 
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Let us demand the energy-momentum tensor to be conserved, i.e., 



t^m — Tf 1 _i_ T@ — T@ — n 

1 u;fj, — 1 i>,tM "r 1 fS/j, 1 v 1 vfj, 1 f3 — u 



(2.21) 



Taking into account that 7T is a function of t only, from Q2.21|) we find 



*«j o -0P o +*ltf-*I^ltf = O, 



c 



(2.22a) 
(2.22b) 



As one can easily verify, the equations ( 2.22[ ) are consistent with those of ( |2.14[ ). 

Let us now deal with the Einstein equations (|2.6| ). In view of (|2.20| ), from (|2.6e| ) one 
obtains 



(2.23) 



whereas, subtraction of ( |2.6b| ) and Q2.6a| ) leads to 



d_ 

dt 



m — n 



(2.24) 



Note that, the two other equations, obtained by subtracting ( |2.6c|) from ( |2.6a| ) and ( |2.6c| ) 
from ( |2.6b| ), respectively, are identical to ( p.24j ), that can be easily verified inserting ( |2.6e|) 
and Q2.14| ) into the equations in question. Finally, summation of ( |2.6a| ), ( |2.6b| ), ( |2.6c| ) and 
( p.6d| ), multiplied by 3, leads to the equation for r, which in view of ( f2.20| ) takes the form 



T - = 2 

T 



in 2 — mn — n 2 



- -[3(MS + VS + QP- 2F) + 2-V 3 }. 
2 c 



;2.25) 



For the right-hand side of the equation ( |2.25| ) are some functions of r only, the solution to 
this equation is well known |15|]. It can be shown that the quantities, related to the spinor 
field are indeed the functions of r. Now, if we assume the metric function c also to be a 
function of r, both ( |2.24| ) and ( |2.25[ ) can be solved explicitly. Thus both the nonlinear spinor 
and the gravitational field equations are solved in general. In what follows, we consider some 
special cases and compare the solutions obtained with those for a BI universe. In can be 
emphasized that the introduction of inhomogeneity both in gravitational (through m and 
n) and spinor (through k) significantly complicates the whole picture. 
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To begin with we consider the linear case setting F(I, J) = 0. As one sees, the structures 
of the spinor field equation ( 2.13Q , as well as the gravitational ones ( |2.24|) and (|2.25|) in this 
case remain unaltered. Thus, the removal of the nonlinearity has little to offer in our 
cause. As it was mentioned earlier, the introduction of inhomogeneity is in the root of all 
these troubles. So, for some break-through we demand the spinor field completely space- 
independent setting k = 0. In this case for the components of the energy-momentum tensor 
immediately we find 



rpO 



mS — F, 



Tl = Tl 



To = 0. 



vs + gp- f, 



(2.26a) 
(2.26b) 
(2.26c) 



In view of ( [2.26| ) from (|2.23|) we obtain 



a m b -n/ c rn-n = ^ 



[2.27) 



but the equations ( |2.13| ), Q2.24j) and ( 2.25|) are still unchanged. Thus we see that even the 
elimination of spinor field inhomogeneity has little to offer. 

Finally, we consider the case, when the spinor field is independent of space i.e., k = 
and the gravitational field is given by a BV space-time with m = n in (|2.4j) . In this case for 
the spinor field we find 



iii + i$u\ — Quz = 0, 
u 2 + i&u 2 - Qu A = 0, 
■u 3 - i&u 3 + Qu\ = 0, 
m 4 — -i$M 4 + Qu 2 = 0. 



(2.28a) 
(2.28b) 
(2.28c) 
(2.28d) 



The bilinear spinor forms in this case satisfy 



S - 2QAI = 0, 
P - 2$A° = 0, 



(2.29a) 
(2.29b) 
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A° + 2$P + 2GS 



A 3 



V° 



V 3 + 2$Q 3 1 ° - 2QQf 



Qo 1 + 20V O 3 



0. 
0. 

0, 

o, 

0. 
0. 



(2.29c) 
(2.29d) 
(2.29e) 
(2.29f) 
(2.29g) 
(2.29h) 



with the relations 



(So) 2 + (P ) 2 + (A c c 



A 3 



1 o 



W, 3 ) 2 + (<2o°) 2 + (Qo 1 ) 2 

with Bi being the constant of integration. 
For the gravitational field we find, 

a = (A0 (1/m) &, 



B 2 , 
B 3 , 
_B 4 , 



(2.30a) 
(2.30b) 
(2.30c) 
(2.30d) 



(2.31) 



with the equations 



and 



d 
dt 



2m 2 



- = -2^- - -[3(MS + VS + QP- 2F)}. 



(2.32) 



(2.33) 



Note that the spinor field equation (|2.28|) completely coincides with those for a BI metric. 
If the nonlinear term in the Lagrangian is given as F = F(I), then the components of the 



spinor field can be given as [12 



Mt) = (Ci/v^)exp [-i / (M - V)dt\ 



(2.34a) 
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Mt) = (CVv^exp [-i f(M- V)dt], (2.34b) 
Mt) = (C 3 /v^)exp [ij(M- V)dt], (2.34c) 
Mt) = (C 4 /Vr)exp [i y*(M - (2.34d) 

with Ci, C2, C3, C4 being the integration constants, such that 

r<2 1 ^t2 ^-y2 ^2 _ s~t 
°1 L/ 2 ~~ °3 ~~ °4 ~~ °0) 

with 

5 = C /r. 

In case, the nonlinear term is given by F = F(J), the components of the spinor field have 
the form 

V>i = -^{Die la + iD 3 e- ia ), (2.35a) 
V r 

^2 = -p^e" + iD 4 e- ia ), (2.35b) 

V> 3 = -^(iDie™ + D 3 e" JCT ), (2.35c) 

ij 4 = -^(iD 2 e ia + D 4 e" !ff ). (2.35d) 
V T 

Here a = intQdt, and the integration constants Di obey 

2 (D? + D\ - D\ - Dl) = D , 
with D to be determined from 

P = D /t. 

Contrary to the BI metric, where the metric functions a, b, c and r are easily deter- 
mined given the concrete form of nonlinearity, in the case considered, the process is much 
complicated with the space inhomogeneity being an active player. Thus we see that the 
introduction of inhomogeneity in the space-time has a far reaching effect in the evolution of 
both spinor and gravitational field. 
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